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Abstract 

Despite the long intertwined evolutionary histories of bees and plants, bee diversity peaks in the xeric areas 
of the eastern and western hemispheres and not the tropics, where plant diversity is greatest. Intensive 
sampling in the northeast Chihuahuan Desert of Mexico and the United States provide the first quan- 
titative estimate of bee species richness where high diversity had been predicted in North America from 
museum records. We find that the density of bee species in a limited area of 16 km” far exceeds any other 
site in the world and amounts to approximately 14% of the bee species described from the United States. 
Long-term studies of bees and other pollinators from areas that are minimally impacted by humans pro- 
vide much-needed baseline data for studies of bees where human impacts are more severe and as climate 


change accelerates. 
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Introduction 


The keystone role of bees in ecosystems has led to greater appreciation and a deeper 
understanding of some previously understudied aspects of their biology. These include 
the contribution of non-managed bee species to crop pollination (Kremen et al. 2002; 
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Winfree et al. 2008; Garibaldi et al. 2013), ecological and life history variables that pre- 
dict vulnerability to environmental change among species (Winfree et al. 2009; Williams 
et al. 2010), and the environmental conditions that favor increases in bee populations 
and their species richness (Cane et al. 2006; Williams et al. 2010; Ollerton et al. 2014). 
Studies that have compared historical records from museums to more recent samples 
have concluded that bees along with other pollinators are in decline due to human activ- 
ity (Biesmeijer et al. 2006; Lebuhn et al. 2012; Burkle et al. 2013; Ollerton et al. 2014; 
Bartomeus et al. 2018; Mathiasson and Rehan 2019). Comparisons of bee faunas at sites 
where land-use histories differ have reached similar conclusions (reviewed in Winfree 
et al. 2011: but see Winfree et al. 2007; Minckley 2014). However, with some notable 
recent exceptions (Herrera 2019; Meiners et al. 2019 and references therein), most of 
this recent work has not extended beyond one or several years and has been focused on 
faunas in or near human-disturbed habitats. Needed are more long-term studies from 
habitats that are little disturbed by humans to provide baseline data to appropriately 
gauge how bee dynamics are influenced where there is anthropogenic disturbance (Win- 
free 2010; Archer et al. 2014; De Palma et al. 2016). Herrera (2019) has gone further to 
emphasize that studies are needed particularly where diversity is highest. 

This study was done to thoroughly document the number of bee species (=species 
richness) in a bee biodiversity hotspot and serves to address the concern of Herrera 
(2019) mentioned above. We use repeat sampling to estimate the total number of 
bee species (observed as well as those not captured) from a defined area so both the 
number of species (species richness, an area free estimate) and number of species from 
a known area (species density) can be compared to data from other studies that are 
focused on elucidating biogeographic patterns of bees. The amount of area sampled 
strongly predicts species richness in all groups (Rozenzweig 1995) and allows species 
density to be evaluated separately from species richness (Gotelli and Colwell 2001). 
Global patterns of species richness remain one of the central questions to evolutionary 
ecology (Colwell and Lees 2000; Lomolino et al. 2017) and anomalous diversity pat- 
terns sometimes provide unusual insights into the mechanisms that generate diversity 
(Kindlmann et al. 2012). Bees are a species-rich group (N = 20,000 species) and have 
very different and counterintuitive biogeographic patterns from their hosts, the flower- 
ing plants. Using species lists and museum records, Michener (1979) showed that the 
warm desert areas of North America and xeric regions around the Mediterranean Sea of 
the eastern hemisphere were unusually rich in bee species. Why bee diversity does not 
follow the temperate to tropic increase in species richness of their floral hosts remains 
poorly understood (Orr et al. 2020). One hypothesis is the persistently wet soils of 
tropical regions favor fungi and other pathogens of bees that nest in the soil (Michener 
1979). This is consistent with the observation that many bee groups that occur in 
the tropics have their nests above-ground in trees and broken stems where moisture 
is reduced (Michener 2007). Bee diversity in the tropics may also be depressed be- 
cause many tropical groups are highly eusocial with colonies that can be large. If floral 
resources are limiting, at least sometimes, social bees may successfully outcompete 
solitary bees. Also possible is that some feature of the ecosystems where warm xeric 
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climates occur favor bee diversification. In parts of the Chihuahuan and Sonoran De- 
serts of North America, and around the Mediterranean Sea in the eastern hemisphere 
rainfall and temperature patterns result in separate spring and summer blooms. Where 
this occurs, bee activity is also biseasonal and many of these bee species are active in 
only one season. The temporally separated but geographically co-occurring bee fauna 
may allow more species to occur in one area than could coexist otherwise (Michener 
1979; Minckley 2008). Finally, topographical complexity and biogeography may also 
contribute to global patterns of bee diversity. Where bee species richness peaks are also 
areas with ongoing mountain building and associated tectonics (Casas-Sainz and de 
Vicente 2009; Gonzalez-Leén et al 2011) that have been shown to favor diversification 
in other groups (Rahbek et al. 2019), and are where biomes converge (Spector 2002). 
Although it is possible that one factor is responsible for most variation in bee diversity, 
it is more likely that some combination of the factors listed above underlie where and 
how many bees occur worldwide. 

We examined the hypothesis that bee richness peaks in the warm deserts of North 
America using intensive, standardized methods in one area of the Chihuahuan Desert 
along the international border between the United States and Mexico. Where this 
study was done had been cattle ranchland since the early 1800's, but has not been 
grazed for the past 20 years in Mexico and the past 50 years in the United States 
(Minckley 2013). This region is a biogeographical crossroads where the subtropical 
dry forest, Great Plains, Chihuahuan and Sonoran Desert intersect (Brown 1994). 
The pollinator fauna is largely native and strongly dominated by bees. The honey bee, 
Apis mellifera, is the sole introduced bee species representing 3.5% of all bee captures. 
The vast majority of collections in this study were within an area of approximately 16 
km? on the floor of the San Bernardino Valley. Repeat sampling using the same pro- 
tocol allowed us to use abundance data to make a quantitative estimate of bee species 
density including species not captured. We then compared estimates of bee species 
density from this 1-year intensive sampling to nine years of collecting using sampling 
by several methods from this same area. Combining sampling methods that vary in 
effort and efficacy violates some assumptions of the models used to estimate the rich- 
ness of observed plus unsampled species (Gotelli and Colwell 2001). However, that 
bee communities, particularly those in deserts, vary dramatically across all spatial and 
temporal scales (Minckley et al. 1999; Williams et al. 2001; Russo et al. 2015; Herrera 
2019) and in response to sampling effort, protocol, and technique (Toler et al. 2005; 
Roulston et al. 2007; Westphal et al. 2008; Droege et al. 2010; Cane et al. 2013; 
Prendergast et al. 2020), suggests that multiple approaches will yield more bee species, 
and more accurate estimates of the number of species. For an unusually well-studied, 
species-rich community of ants in Costa Rica, Longino et al. (2002) showed that the 
actual species richness was greatly underestimated with datasets from single sampling 
protocols; and only by combining different datasets was it possible to obtain an accu- 
rate estimate. In this study, bee species density estimates based on repeat samples with 
one protocol is compared to estimates based on a larger dataset from the same area that 
combined sampling protocols from more sites over nine years. 
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Our prediction was that intensive sampling using one technique (pan traps) over 
one year would underestimate the observed species richness obtained over more years 
using multiple sampling techniques from the same area. However, we also predicted 
this intensive pan trapping would provide a reasonable estimate of species richness 
using statistical methods that extrapolate from species observed to include species not 
captured but present (Gotelli and Colwell 2001). 

We also gauged how well our samples from the valley floor estimate the number of 
species in a broader area, by compiling a dataset from collections made in two nearby 
areas and making an assumption that this list of species from all collections combined 
was the best estimate for bee richness in this region. Extrapolation to the expected spe- 
cies richness from the multiyear samples made on the valley floor closely approximated 
the observed species richness of the larger regional dataset. In comparison to 21 other 
long-term studies of bees, the density of bee species in this part of the Chihuahuan 
Desert is greater than that reported from any other area. 


Methods 


The study was done in the San Bernardino Valley in the northeastern Chihuahuan De- 
sert of western North America near the confluence of Arizona and New Mexico, USA 
and Chihuahua and Sonora, Mexico. The San Bernardino Valley runs north-south 
across the Mexico-United States border in northeastern Sonora, Mexico and southeast- 
ern Arizona, USA. Elevation is approximately 1070 m and climate is xeric temperate 
with an annual average precipitation of 360 mm/ year. Approximately 50% of the an- 
nual precipitation occurs in July and August and the two driest months are April and 
May based on long-term weather records 30 kilometers west of the San Bernardino 
Valley at Douglas, Arizona. Minckley (2008, 2013) provides a detailed description of 
the climate, area and vegetation. 

The main study area contains a number of habitats ranging from desert scrub to 
permanent springs that are the headwaters of the San Bernardino River. In 2000, seven 
to eleven permanently marked 1-ha sites were established in five different vegetation 
types (Fig. 1). Two vegetation types were in proximity to surface (riparian) or close 
subsurface (mesquite forest) water and were more mesic than other vegetation types. 
Two vegetation types were upland and were much drier (grassland, desert scrub) with 
cacti and shrubs characteristic of Chihuahuan Desert, such as Acacia constricta Benth. 
(cat claw), Larrea tridentata (DC.) Coville (creosote bush), and Flourensia cernua DC. 
(tarbush). One habitat was abandoned agricultural fields. All sites were within a 16 
km? area and ranged from 1100 to 1220 m in elevation. 

Collections of bees were made at all of the permanent sites by pan trap in 2001 
and at approximately half of these sites from 2002 to 2008. Pan trap sampling was 
done using three yellow, three white and three blue 177 ml. pan traps placed out in 
alternating colors and spaced 3-4 m apart along a transect following the protocol of 


LeBuhn et al. (2007). 
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Figure |. Map of the study area in the San Bernardino Valley (modified from Minckley 2014). Dots 
are the permanent sites where 2001 pan trap collections were made. Letters next to each permanent site 
indicate the habitat; F = field, G = grassland, M = mesquite, R = riparian, and S = scrub. The core area 
is delimited by the rectangle spanning the Mexico — United States border. The peripheral and limestone 
collection sites are not shown but the direction and distance from the center of the core area are indicated 
on the right side of the map. 


All bee specimens were identified to species and sex by taxonomic experts and pub- 
lished keys (Suppl. material 1: Table $1). Most were identified to recognized species, 
except for 73 species that were recognized as morphospecies. These morphospecies are 
1.43% of all individuals, and most are members of Lasioglossum (Dialictus), Dufourea 
and Megachile. All species and their abundance are listed in Table 1 and are held in the 
Minckley collection at the University of Rochester. Data are available from the Global 
Biodiversity Information Facility (GBIF) under the RL Minckley Insect and Plant 
Collection (https://www.gbif.org/dataset/84b69033-f9de-40bc-8d5f-8ee41987d6cc). 


Species density 


Bee species density in the San Bernardino Valley was first estimated from the bee col- 
lections captured in pan traps at permanent sites in 2001 (Fig. 1). These samples were 
made at all 45 permanent sites on 10 to 14-day intervals from mid-April to mid- 
September, except for three weeks in June, the driest time of the year (Tables 2, 3). 

A second estimate of bee species density was done from all samples in the same 
area circumscribed by the permanent sites described above (hereafter referred to as the 
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Table |. List of bee species and their abundance found in the San Bernardino Valley, Sonora/Arizona. 
Columns refer to numbers of individuals from the four sets of samples: A is the 2001 pan trap samples 
made at permanent sites; B is the core area samples taken from 2000 to 2008 from the area circumscribed 
by the permanent sites less those samples in A; C is the samples taken from 2000 to 2008 in the area 
peripheral to the core area; D is samples taken from 2009-2014 from a limestone outcrop area of desert 
approximately 12 km north northeast of the valley floor area. The taxonomy follows that of Michener 
(2007) except for recent changes by Bossert et al. (2020). 


A B C D 
Andrenidae 
Andreninae 
Ancylandrena rozeni Zavortink, 1994 0 0 2 0 
Andrena (Belandrena) sphaeralceae Linsley, 1939 0 0 Di 0 
A, (Diandrena) olivacea Viereck, 1917 0 1 0 0 
A. (Euandrena) auricoma Smith, 1879 0 16 4 0 
A. (Holandrena) cressonii Robertson, 1891 0 0 0 1 
A. (Leucandrena) monilicornis Cockerell, 1896 0 1 3 1 
A, (Melandrena) cerasifolii Cockerell, 1896 0 2 0 0 
A. (Micrandrena) piperi Viereck, 1904 2 28 1 0 
A, (Plastandrena) fracta Casad & Cockerell, 1896 4 111 34 19 
A. (Plastandrena) prunorum Casad, 1896 0 1 0 59 
A. (Rhaphandrena) prima Casad, 1896 0 3 1 Z 
A. (Scaphandrena) capricornis Casad & Cockerell, 1896 0 6 2 13 
A. (Scaphandrena) primulifrons Casad, 1896 3 83 12 0 
A. (Thysandrena) w-scripta Viereck, 1904 0 0 0 51 
A. (Trachandrena) semipunctata Cockerell, 1902 0 9 2 0 
A, (Tylandrena) jessicae Cockerell, 1896 0 1 0 0 
A, sp. 4 0 0 0 3 
A. sp. 10 3 8 1 0 
A. sp. 11 0 1 2 0 
A. sp. 19 0 1 12 0 
A. sp. 23 0 1 3 0 
A. sp. 27 0 10 1 0 
A. sp. 30 0 6 0 0 
A. sp. 36 1 5 0 0 
Oxaeinae 
Protoxea gloriosa (Fox, 1893) 3 29 0 0 
Panurginae 
Calliopsini 

Calliopsis sonorana (Timberlake, 1969) 10 1 0 0 
C. (Calliopsima) chlorops Cockerell, 1899 0 1 0 0 
C. (Calliopsima) coloratipes Cockerell, 1898 0 4 0 0 
C. (Calliopsima) crypta Shinn, 1965 0 1 0 0 
C. (Calliopsima) pectidis Shinn, 1965 0 1 0 0 
C. (Calliopsima) rozeni Shinn, 1965 41 46 3 0 
C. (Calliopsis) empelia Shinn, 1967 4 0 0 
C. (Calliopsis) sonora Shinn, 1967 3 > 0 0 
C. (Calliopsis) squamifera Timberlake, 1947 1 13 0 0 
C. (Hypomacrotera) callops (Cockerell & Porter, 1899) 207 482 27 0 
C. (Hypomacrotera) persimilis (Cockerell, 1899) 597 1071 11 0 
C. (Hypomacrotera) subalpina subalpina Cockerell, 1894 23 257 32 0 
C. (Nomadopsis) australior Cockerell, 1897 8 70 0 0 
C. (Nomadopsis) callosa Vimberlake, 1952 18 8 0 0 
C. (Nomadopsis) helianthi Swenk & Cockerell, 1907 36 176 3 0 
C. (Nomadopsis) macswaini (Rozen, 1958) 0 6 0 0 
C. (Nomadopsis) nigromaculata Timberlake, 1952 3 4 0 0 
C. (Nomadopsis) obscurella Cresson, 1879 0 1 4 0 


Bee diversity hotspot 323 


A B C D 
C. (Nomadopsis) puellae (Cockerell, 1933) 1019 752 61 2 
C. (Perissander) anomoptera Michener, 1942 4 16 1 0 
C. (Perissander) gilva Shinn, 1967 ik 19 1 0 
C. (Perissander) limbus Shinn, 1967 0 1 0 0 
C. (Perissander) rogeri Shinn, 1967 2 2 1 0 
C. (Verbenapis) verbenae Cockerell & Porter, 1899 4 2 0 0 

Perditini 
Macrotera (Cockerellula) parkeri (Timberlake, 1980) 0 0 0 216 
M. (Cockerellula) rubida (Yimberlake, 1968) 1 0 0 0 
M. (Cockerellula) solitaria (Cockerell, 1897) 2 1 0 0 
M. (Macrotera) texana Cresson, 1878 0 1 0 0 
M. (Macroterella) mellea Timberlake, 1954 2 16 0 0 
M. (Macroteropsis) latior (Cockerell, 1896) 78 123 10 0 
M. (Macroteropsis) portalis (Timberlake, 1959) 6 41 0 0 
Perdita (Cockerellia) albipennis Cresson, 1868 0 1 0 0 
P (Cockerellia) coreopsidis collaris Cockerell, 1916 0 2 0 0 
P (Cockerellia) verbesinae Cockerell, 1896 0 0 0 
P (Epimacrotera) biguttata Timberlake, 1962 60 10 0 0 
P (Epimacrotera) diversa Timberlake, 1954 315 313 11 0 
P. (Glossoperdita) hurdi Timberlake, 1956 0 28 0 
P (Heteroperdita) rhodogastra Timberlake, 1954 35 68 0 
P. (Hexaperdita) asteris Cockerell, 1986 3 0 0 0 
P. (Hexaperdita) callicerata Cockerell, 1896 238 360 30 0 
P. (Hexaperdita) ignota Cockerell, 1896 3 1 0 0 
P (Pentaperdita) albovittata Cockerell, 1895 59 180 1 0 
P. (Pentaperdita) amoena Timberlake, 1956 173 67 5 0 
P (Pentaperdita) bradleyana Timberlake, 1954 0 4 0 0 
P (Pentaperdita) melanochlora Cockerell, 1922 2 0 0 0 
P (Perdita) affinis Cresson, 1878 1 0 0 0 
P (Perdita) ashmeadi Cockerell, 1899 1 196 9 0 
P (Perdita) chamaesarachae Cockerell, 1896 3 31 16 0 
P (Perdita) dasylirii Cockerell, 1907 0 0 46 0 
P (Perdita) dificilis Timberlake, 1964 0 3 0 0 
P (Perdita) exclamans Cockerell, 1895 8 301 30 0 
P (Perdita) florisantella Cockerell, 1906 47 148 38 0 
P (Perdita) lenis Timberlake, 1958 3 232 46 0 
P (Perdita) luciae Cockerell, 1899 0 50 0 0 
P (Perdita) mimosae efferta Timberlake, 1964 1 0 0 0 
P (Perdita) munita Timberlake, 1964 205 337 20 0 
P (Perdita) pectidis Cockerell, 1896 6 13 0 0 
P (Perdita) punctifera Cockerell, 1914 0 161 0 0 
P (Perdita) punctosignata Cockerell, 1895 5 112 1 0 
P (Perdita) semicaerulea Cockerell, 1896 1 467 21 0 
P (Perdita) semicrocea Cockerell, 1895 5 1 0 0 
P (Perdita) sexmaculata Cockerell, 1895 31 76 81 0 
P (Perdita) stathamae Timberlake, 1964 0 43 8 0 
P (Perdita) triangulifera Timberlake, 1964 2, 59 16 0 
P (Perditella) cladothricis Cockerell, 1896 15 23 0 0 
P (Perditella) larreae Cockerell, 1896 18 449 1 0 
P (Perditella) minima Cockerell, 1923 8 9 0 0 
P. (Pygoperdita) malacothricis Timberlake, 1956 887 454 26 3 
Perdita sp. 2 0 29 0 0 
Protandrenini 

Protandrena (Heterosarus) nanula (Timberlake, 1964) 4 18 0 0 
P (Heterosarus) townsendi (Cockerell, 1897) 2 4 0 0 
P (Protandrena) bancrofti Dunning, 1897 1 3 0 0 
Pseudoanurgus fraterculus fraterculus Cockerell, 1896 1 0 0 0 
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Apidae 
Anthophorinae 


Anthophorini 


Anthophora (Anthophoroides) californica Cresson, 1869 


A. (Anthophoroides) vallorum (Cockerell, 1896) 
A, (Lophanthophora) affabilis Cresson, 1878 

A, (Lophanthophora) ursina Cresson, 1869 

A. (Melea) bomboides Kirby, 1838 

A. (Micranthophora) aff. estebana 

A, (Micranthophora) curta Provancher, 1895 

A, (Micranthophora) exigua Cresson, 1878 

A, (Micranthophora) pachyodonta Cockerell, 1923 
A. (Micranthophora) aff. estebana 

A, (Heliophila) petrophila Cockerell, 1905 

A, (Mystacanthophora) montana Cresson, 1869 
A. (Mystacanthophora) urbana Cresson, 1878 
A, (Paramegilla) centriformis Cresson, 1879 

A, (Paramegilla) fulvicauda Timberlake, 1937 
A, (Pyganthophora) lesquerellae Cresson, 1878 
A, (Pyganthophora) vannigera Vimberlake, 1951 


Apinae 


Apini 


Apis (Apis) mellifera Linnaeus, 1758 


Bombini 


Bombus (Fervidobombus) pensylvanicus sonorus Say, 1837 


Centridini 


Centris (Paracentris) angustifrons Snelling, 1966 
C. (Paracentris) atripes Mocsary, 1899 

C. (Paracentris) caesalpiniae Cockerell, 1897 

C. (Paracentris) hoffmanseggiae Cockerell, 1897 
C. (Paracentris) rhodopus Cockerell, 1897 

C. (Paracentris) apache Vivallo, 2020 

C. (Paracentris) sp. 


Eucerinae 


Ancyloscelidini 


Ancyloscelis melanostomus Cockerell, 1923 


Emphorini 


Diadasia (Coquillettapis) afflictula Cockerell, 1910 
D. (Coquillettapis) australis (Cresson, 1878) 

D. (Coquillettapis) diminuta (Cresson, 1878) 

D. (Coquillettapis) lutzi Cockerell, 1924 

D. (Coquillettapis) martialis Timberlake, 1940 

D. (Coquillettapis) megamorpha Cockerell, 1898 
D. (Coquillettapis) opuntiae Cockerell, 1901 

D. (Coquillettapis) rinconis Cockerell, 1897 

D. (Coquillettapis) sphaeralcearum Cockerell, 1897 
D. (Dasiapis) ochracea (Cockerell, 1903) 

D. (Dasiapis) tropicalis (Cockerell, 1918) 
Ptilothrix nr. sumichrasti (Cresson, 1878) 


Eucerini 


Eucera (Synhalonia) actuosa (Cresson, 1878) 

E. (Synhalonia) albescens (Timberlake, 1969) 

E. (Synhatonia) aragalli (Cockerell, 1904) 

E. (Synhalonia) conformis (Timberlake, 1969) 

E. (Synhalonia) edwardsi (Cresson, 1878) 

E. (Synhalonia) lepida (Cresson, 1878) 

E. (Synhalonia) mohavensis (Timberlake, 1969) 

E. (Synhalonia) nr. quadricincta (Timberlake, 1969) 
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E. (Synhalonia) phaceliae (Cockerell, 1911) 

E. (Synhalonia) quadricincta (Timberlake, 1969) 

E. (Synhatonia) sp. 1 

E. (Synhatonia) sp. 2 

Martinapis (Martinapis) luteicornis (Cockerell, 1896) 

Melissodes (Eumelissodes) agilis Cresson, 1878 

M. (Eumelissodes) limbus LaBerge, 1961 

M. (Eumelissodes) lutulentus LaBerge, 1961 

M. (Eumelissodes) subagilis Cockerell, 1905 

M. (Eumelissodes) tristis Cockerell, 1894 

M. (Eumelissodes) verbesinarum Cockerell, 1905 

M. (Melissodes) communis Cresson, 1878 

M. (Melissodes) comptoides Robertson, 1878 

M. (Melissodes) gilensis Cockerell, 1896 

M. (Melissodes) paroselae Cockerell, 1905 

M. (Tachymelissodes) opuntiellus Cockerell, 1911 

M. (Tachymelissodes) sonorensis LaBerge, 1963 

M. sp. 22 

Peponapis (Peponapis) pruinosa (Say, 1837) 

P. (Xeropeponapis) timberlakei Hurd & Linsley, 1964 

Svastra (Epimelissodes) helianthelli (Cockerell, 1905) 

S. (Epimelissodes) machaerantherae (Cockerell, 1904) 

S. (Epimelissodes) obliqua expurgata (Cockerell, 1925) 

S. (Epimelissodes) sabinensis (Cockerell, 1924) 

S. (Epimelissodes) sila (LaBerge, 1956) 

S. (Epimelissodes) sp. 4 

Syntrichalonia exquisita (Cresson, 1878) 

Tetraloniella (Tetraloniella) eriocarpi (Cockerell, 1898) 

T. (Tetraloniella) imitatrix (Cockerell & Porter, 1899) 

T. (Tetraloniella) lippiae (Cockerell, 1904) 

T. (Tetraloniella) sphaeralceae LaBerge, 2001 

T. (Tetraloniella) vandyckei LaBerge, 2001 

Xenoglossa (Eoxenoglossa) strenua (Cresson, 1878) 

X. (Xenoglossa) angustior Cockerell, 1900 

X. (Xenoglossa) patricia Cockerell, 1896 
Exomalopsini 

Anthophorula (Anthophorisca) asteris (Mitchell, 1962) 

A. (Anthophorisca) consobrina (Timberlake, 1980) 

A, (Anthophorisca) exilis (Timberlake, 1980) 

A, (Anthophorisca) ignota (Timberlake, 1980) 

A, (Anthophorisca) parva (Timberlake, 1980) 

A, (Anthophorisca) pygmaea (Cresson, 1872) 

A, (Anthophorisca) sonorensis (Timberlake, 1980) 

A. (Anthophorula) albata (Timberlake, 1947) 

A, (Anthophorula) compactula (Cockerell, 1897) 

A, (Anthophorula) completa (Cockerell, 1935) 

A, (Anthophorula) crenulata (Timberlake, 1980) 

A. (Anthophorula) euphorbiae (Timberlake, 1947) 

A, (Anthophorula) gutierreziae (Timberlake, 1947) 

A. (Anthophorula) palmarum (Timberlake, 1947) 

A, (Anthophorula) rozeni (Timberlake, 1980) 

A, (Anthophorula) rufiventris (Timberlake, 1947) 

A, (Anthophorula) scapalis (Timberlake, 1980) 

A, (Anthophorula) torticornis (Timberlake, 1980) 

A, (Anthophorula) tricinctula (Timberlake, 1980) 

A, (Anthophorula) varleyi (Timberlake, 1980) 

Exomatopsis (Stilbomalopsis) dimidiata Timberlake, 1980 

E. (Stilomatopsis) solani Cockerell, 1896 
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E. (Stilomalopsis) solidaginis Cockerell, 1898 


Nomadinae 


Ammobatini 


Oreopasites (Oreopasites) arizonica Linsley, 1941 
O. (Perditopasites) favreauae Rozen, 1992 
O. (Perditopasites) linsleyi Rozen, 1992 


Ammobatoidini 


Holcopasites apacheorum (Hurd & Linsley, 1972) 
LH. illinoiensis minimus (Linsley, 1943) 

H. insoletus (Linsley, 1942) 

H. stevensi Crawford, 1915 

H. tegularis Hurd & Linsley, 1972 


Brachynomadini 


Brachynomada (Melanomada) margaretae (Rozen, 1994) 


Paranomada velutina Linsley, 1939 
Triopasites penniger (Cockerell, 1894) 


Epeolini 


Epeolus chamaesarachae Onuferko, 2018 
E. compactus Cresson, 1878 

E. mesillae (Cockerell, 1895) 

Epeolus sp. 1 

Epeolus sp. 3 

Triepeolus circumculus Rightmyer, 2008 
T. concavus (Cresson, 1878) 

T. grandis (Friese, 1917) 

T. helianthi (Robertson, 1897) 

T. kathrynae Rozen, 1989 

T. loomisorum Rozen, 1989 

T. lunatus (Say, 1824) 

T. norae Cockerell, 1907 

T. penicilliferus (Brues, 1903) 

T. remigatus (Fabricius, 1804) 
Triepeolus sp. 4 

T. townsendi Cockerell, 1907 

T. verbesinae Cockerell, 1897 


Ericrocidini 


Ericrocis lata (Cresson, 1878) 
E. pintada Snelling & Zavortink, 1984 


Melectini 


Melecta (Melecta) bohartorum Linsley, 1939 


Neolarrini 


Neopasites (Micropasites) cressoni Crawford, 1916 
Neolarra (Neolarra) batrae Shanks, 1977 

N. (Neolarra) californica Michener, 1939 

N. (Neolarra) cockerelli (Crawford, 1916) 

N. (Phileremulus) rozeni Shanks, 1977 

N. (Phileremulus) vigilans (Cockerell, 1895) 
Townsendiella pulchra Crawford, 1916 


Nomadini 


Nomada cf. gutierreziae Cockerell, 1896 
Nomada cf. vegana Cockerell, 1903 
Nomada sp. 1 

Nomada sp. la 

Nomada sp. 2 

Nomada sp. 3 

Nomada sp. 4 

Nomada sp. 5 

Nomada sp. 6 
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A B C D 
Nomada sp. 7 0 3 1 1 
Nomada sp. 8 0 1 0 0 
Nomada sp. 9 0 0 0 1 
Nomada sp. 10 0 0 0 3 
Nomada sp. 11 0 3 1 0 
Xylocopinae 
Ceratinini 
Ceratina (Zadontomerus) apacheorum Daly, 1973 2324 2849 11 2 
C. (Zadontomerus) arizonensis Cockerell, 1898 8 188 
C. (Zadontomerus) melanoptera Cockerell, 1924 34 5 3 15 
C. (Zadontomerus) nanula Cockerell, 1897 0 1 0 0 
Xylocopini 
Xylocopa (Xylocopoides) californica arizonensis Cresson, 1879 0 8 oD: 2 
X. (Neoxylocopa) sonorina Smith, 1874 0 2 0 0 
Colletidae 
Colletinae 
Colletini 
Colletes bryanti Timberlake, 1951 0 2 0 0 
C. clypeonitens Swenk, 1906 2 3 0 0 
C. eulophi Robertson, 1891 0 2 0 0 
C. intermixtus Swenk, 1905 0 1 0 0 
C. kincaidti Cockerell, 1898 0 1 0 0 
C. louisae Cockerell, 1897 0 12 1 0 
C. lutzi Timberlake, 1943 0 1 0 0 
C. perileucus Cockerell, 1924 0 18 0 0 
C. prosopidis Cockerell, 1897 0 2 0 0 
C. simulans simulans Cresson, 1868 0 2 0 0 
C. salicicola Cockerell, 1897 1 141 17 0 
C. scopiventer Swenk, 1908 0 57 17 0 
C. sphaeralceae Timberlake, 1951 0 1 1 0 
C. wootoni Cockerell, 1897 0 74 18 0 
Diphaglossinae 
Caupolicanini 
Caupolicana (Caupolicana) yarrowi (Cresson, 1875) 0 3 0 0 
C. (Zikanapis) elegans Timberlake, 1965 0 0 1 0 
Hylaeinae 
Hylaeus (Hylaeus) mesillae (Cockerell, 1896) 18 469 5 0 
H. (Paraprosopis) asininus (Cockerell & Casad, 1895) 2 90 11 0 
H. (Prosopis) episcopalis (Cockerell, 1896) 0 340 0 
H. (Prosopis) insolitus Snelling, 1966 0 1 2 0 
Halictidae 
Halictinae 
Augochlorini 
Augochlorella (Augochlorella) neglectula (Cockerell, 1897) 4 21 27 6 
A, (Augochlorella) pomoniella (Cockerell, 1915) 2 4 2 4 
Halictini 
Agapostemon (Agapostemon) angelicus Cockerell, 1924 200 684 67 ol 
A, (Agapostemon) melliventris Cresson, 1874 32 75 0 11 
A, (Agapostemon) obliquus (Provancher, 1888) 402 3 11 78 
A, (Agapostemon) texanus Cresson, 1872 0 0 0 
A, (Agapostemon) tyleri Cockerell, 1917 18 153 7 12 
A, (Notagapostemon) nasutus Smith, 1853 0 1 0 0 
Augochloropsis metallica (Fabricius, 1793) 0 12 2 0 
Halictus (Odontalictus) ligatus Say, 1837 31 578 31 0 
H. (Seladonia) tripartitus Cockerell, 1895 1296 1042 63 D) 
Lasioglossum (Dialictus) arcanum (Sandhouse, 1924) 0 1 0 
L. (Dialictus) clematisellum (Cockerell, 1904) 8 252 3 0 


L. (Dialictus) comulum Michener, 1951 401 768 81 675 
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L. (Dialictus) eophilus (Ellis, 1914) 

L. (Dialictus) hudsoniellum (Cockerell, 1919) 
L. (Dialictus) impavidum (Sandhouse, 1924) 
L. (Dialictus) lilianae 

L. (Dialictus) mesillense 

L. (Dialictus) microlepoides (Ellis, 1914) 

L. (Dialictus) minckleyi Gardner & Gibbs, 2020 
L. (Dialictus) obnubilum (Sandhouse, 1924) 
L. (Dialictus) occidentale (Crawford, 1902) 

L. (Dialictus) perparvum (Ellis, 1914) 

L. (Dialictus) pruinosum (Robertson, 1892) 

L. (Dialictus) semibrunneum (Cockerell, 1895) 
L. (Dialictus) semicaeruleum (Cockerell, 1895) 
L. (Dialictus) aff. macroprosopum Gibbs, 2010 
L. (Dialictus) aff. occidentale (Crawford, 1902) 
L. (Dialictus) aff. perparvum (Ellis, 1914) 


L. (Dialictus) cf. albuquerquense (Michener, 1937) 


L. (Dialictus) cf. brunneri (Crawford, 1902) 
L. (Dialictus) cf. comulum (Michener, 1937) 
L. (Dialictus) cf. impavidum (Sandhouse, 1924) 
L. (Dialictus) cf. lionotum (Sandhouse, 1923) 
L. (Dialictus) cf. nevadense (Crawford, 1907) 
L. (Dialictus) cf. obnubilum (Sandhouse, 1924) 
L. (Dialictus) cf. perdifficile (Cockerell, 1895) 
L. (Dialictus) sp. 1 

L. (Evylaeus) amicum (Cockerell, 1897) 

L. (Evylaeus) pectoraloides (Cockerell, 1895) 

L. (Evylaeus) sp. 1 

L. (Evylaeus) sp. 4 

L. (Hemihalictus) angustius (Cockerell, 1897) 
L. (Lasioglossum) acarophilum McGinley, 1986 
L. (Lasioglossum) desertum (Smith, 1879) 

L. (Lasioglossum) morrilli (Cockerell, 1919) 

L. (Lasioglossum) jubatum (Vachal, 1904) 

L. (Lasioglossum) sisymbrii (Cockerell, 1897) 
Sphecodes fortior Cockerell, 1898 

S. minor Robertson, 1898 

S. aff. mandibularis Cresson, 1872 

Sphecodes sp. 1 


Nomiinae 


Dieunomia (Epinomia) boharti (Cross, 1958) 
D. (Epinomia) nevadensis (Cresson, 1874) 
Nomia (Acunomia) angustitibialis Ribble, 1965 
N. (Acunomia) foxii Dalla Torre, 1896 

N. (Acunomia) tetrazonata (Cresson, 1874) 


Rhophitinae 


Conanthalictus (Phaceliapis) sp. 1 
Dufourea matacothricis Timberlake, 1939 
D. mulleri (Cockerell, 1898) 

D. pulchricornis (Cockerell, 1916) 
Dufourea sp. 1 

Dufourea sp. 2 
Dufourea sp. 3 
Dufourea sp. 4 
Dufourea sp. 5 
Dufourea sp. 6 


— 


Dufourea sp. 10 
Dufourea sp. 11 
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A B C D 
Dufourea sp. 12 0 4 Zz 0 
Dufourea sp. 13 0 2 1 0 
Protodufourea eickworti Bohart & Griswold, 1997 5 2 0 0 
Sphecodosoma (Sphecodosoma) dicksoni (Timberlake, 1961) 2 19 0 0 
S. (Sphecodosoma) pratti Crawford, 1907 0 2 0 0 
Megachilidae 
Megachilinae 

Anthidiini 
Anthidium (Anthidium) cochimi Snelling, 1992 6 10 1 0 
A. (Anthidium) cockerelli Schwarz, 1928 8 106 4 0 
A. (Anthidium) jocosum Cresson, 1878 471 390 11 7 
A, (Anthidium) labergei Gonzalez & Griswold, 2013 1 13 6 0 
A, (Anthidium) maculifrons Smith, 1854 3 3 0 0 
A. (Anthidium) maculosum Cresson, 1878 6 13 1 0 
A, (Anthidium) palmarum Cockerell, 1904 34 50 2 1 
A, (Anthidium) paroselae Cockerell, 1898 4 0 1 0 
A, (Anthidium) quetzalcoatli Schwarz, 1933 0 1 0 0 
A. (Anthidium) utahense Swenk, 1914 3 13 2 0 
A. (Anthidium) sp. 4 3 1 1 0 
Dianthidium (Dianthidium) curvatum (Smith, 1854) 0 1 1 0 
D. (Dianthidium) discors Timberlake, 1948 0 3 0 1 
D. (Dianthidium) heterulkei Schwarz, 1940 0 2 1 0 
D. (Dianthidium) implicatum Timberlake, 1948 1 3 0 0 
D. (Dianthidium) parkeri Grigarick & Stange, 1964 0 4 0 0 
D. (Dianthidium) parvum (Cresson, 1878) 6 16 0 0 
D. (Dianthidium) platyurum Cockerell, 1923 0 3 0 0 
D. (Dianthidium) pudicum (Cresson, 1879) 3 1 2 0 
D. (Dianthidium) ulkei (Cresson, 1878) 0 2 0 0 
Stelis (Dolichostelis) perpulchra Crawford, 1916 0 136 4 0 
S. (Stelis) elongativentris Parker, 1987 105 599 6 12 
Stelis sp. 2 0 2 0 0 
Stelis sp. 3 0 2 0 0 
Stelis sp. 4 1 31 1 0 
Stelis sp. 5 2 2 1 0 
Stelis sp. 6 0 16 0 0 
Stelis sp. 7 0 4 1 0 
Stelis sp. 8 0 3 0 0 
Stelis sp. 9 1 2 0 0 
Stelis sp. 10 0 1 0 0 
Trachusa (Heteranthidium) larreae (Cockerell, 1897) 3 137 6 0 

Dioxyini 
Dioxys pomonae Cockerell, 1910 0 0 0 1 
D. producta subrubra (Cresson, 1879) 0 3 0 

Lithurgini 
Lithurgopsis (Lithurgopsis) apicalis (Cresson, 1875) 11 39 24 0 
L. (Lithurgopsis) echinocacti Cockerell, 1898 4 19 0 0 
L. (Lithurgopsis) planifrons (Friese, 1908) 0 3 0 0 

Megachiliini 
Coelioxys (Boreocoelioxys) novomexicanus Say, 1824 2 4 1 0 
Coelioxys (Boreocoelioxys) octodentatus Say, 1824 0 24 1 0 
C. (Coelioxys) hirsutissimus Cockerell, 1912 1 5 3 0 
C. (Syncoelioxys) apacheorum Cockerell, 1900 0 2 0 0 
C. (Syncoelioxys) hunteri Crawford, 1914 0 1 0 0 
C. (Syncoelioxys) texanus Cresson, 1872 0 22 1 0 
C. (Xerocoelioxys) edita Cresson, 1872 1 6 0 0 
Megachile (Argyropile) parallela Smith, 1853 0 2 0 0 
M. (Chelostomoides) adelphodonta Cockerell, 1924 0 5 3 0 
M. (Chelostomoides) chilopsidis Cockerell, 1900 0 31 0 0 
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M. (Chelostomoides) discorhina Cockerell, 1924 
M. (Chelostomoides) exilis Cresson, 1872 

M. (Chelostomoides) lobatifrons Cockerell, 1924 
M. (Chelostomoides) manni Mitchell, 1934 

M. (Chelostomoides) odontostoma Cockerell, 1924 
M. (Chelostomoides) prosopidis Cockerell, 1900 
M. (Chelostomoides) reflexa (Snelling, 1990) 
M. (Chelostomoides) spinotulata Mitchell, 1934 
M. (Leptorachis) petulans Cresson, 1878 

M. (Litomegachile) brevis Say, 1837 

M. (Litomegachile) coquilletti Cockerell, 1915 
M. (Litomegachile) gentilis Cresson, 1872 

M. (Litomegachile) lippiae Cockerell, 1900 

M. (Litomegachile) mendica Cresson, 1878 

M. (Litomegachile) pankus Bzdyk, 2012 

M. (Litomegachile) texana Cresson, 1878 

M. (Megachile) sp 1 

M. (Megachiloides) alata Mitchell, 1934 

M. (Megachiloides) anograe Mitchell, 1938 

M. (Megachiloides) bradleyi Mitchell, 1934 

M. (Megachiloides) fucata Mitchell, 1934 

M. (Megachiloides) integra Cresson, 1878 

M. (Megachiloides) legalis Cresson, 1879 

M. (Megachiloides) maurata Mitchell, 1936 
M. (Megachiloides) sublaurita Mitchell, 1927 
M. (Megachiloides) xerophila Cockerell, 1933 
M. (Megachiloides) sp. 1 

M. (Megachiloides) sp. 2 

M. (Megachiloides) sp. 3 

M. (Pseudocentron) sidalceae Cockerell, 1897 
M. (Sayapis) inimica sayi (Cresson, 1878) 

M. (Sayapis) newberryae Cockerell, 1900 

M. (Sayapis) policaris Say, 1831 

M. (Sayapis) pugnata Say, 1837 


Osmiini 


Ashmeadiella (Arogochila) breviceps Michener, 1939 
A, (Arogochila) cazieri Michener, 1939 

A. (Arogochila) clypeodentata Michener, 1936 
A, (Arogochila) sp. 1 

A. (Ashmeadiella) bigeloviae (Cockerell, 1897) 
A. (Ashmeadiella) bucconis (Say, 1837) 

A. (Ashmeadiella) cactorum (Cockerell, 1897) 
A. (Ashmeadiella) foveata Michener, 1939 

A. (Ashmeadiella) gillettei Titus, 1904 

A. (Ashmeadiella) leucozona Cockerell, 1924 
A. (Ashmeadiella) maxima Michener, 1936 
A. (Ashmeadiella) meliloti (Cockerell, 1897) 
A. (Ashmeadiella) occipitalis Michener, 1939 
A, (Ashmeadiella) opuntiae (Cockerell, 1897) 
A. (Ashmeadiella) prosopidis (Cockerell, 1897) 
A. (Ashmeadiella) sonora Michener, 1939 

A. (Ashmeadiella) vandykiella Michener, 1949 
A. (Ashmeadiella) sp. 2 

A. (Chilosmia) rhodognatha Cockerell, 1924 
A. (lsosmia) hurdiana (Michener, 1954) 

A. (lsosmia) rubrella (Michener, 1943) 
Atoposmia (Atoposmia) anthodyta (Michener, 1943) 
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A B C D 
A. (Atoposmia) arizonensis (Michener, 1954) 2 2 1 0 
A. (Eremosmia) beameri (Michener, 1951) 1 4 1 0 
A. (Evemosmia) hypostomalis (Michener, 1954) 0 1 1 1 
Atoposmia sp. 1 2 4 1 0 
A. (Hexosmia) copelandica arefacta (Cockerell, 1935) 4 19 48 0 
Heriades (Neotrypetes) crucifera Cockerell, 1897 0 4 0 0 
H. (Neotrypetes) micheneri Timberlake, 1947 0 8 0 0 
H. (Neotrypetes) microthalma Michener, 1954 0 5 0 0 
H. (Neotrypetes) texana Michener, 1938 0 1 0 0 
H. (Neotrypetes) timberlakei Michener, 1938 0 2 0 0 
Hoplitis (Alcidamea) grinnelli (Cockerell, 1910) 0 8 4 0 
H. (Alcidamea) biscutellae (Cockerell, 1897) 18 241 13 4 
H. (Alcidamea) producta interior Michener, 1947 0 1 2 0 
H. (Proteriades) remotula (Cockerell, 1910) 0 1 0 0 
Osmia (Diceratosmia) subfasciata Cresson, 1872 291 539 68 4 
O. (Melanosmia) cerasi Cockerell, 1897 0 0 1 0 
O. (Melanosmia) clarescens Cockerell, 1911 2 69 8 0 
O. (Melanosmia) liogastra Cockerell, 1933 3 46 16 34 
O. (Melanosmia) marginata Michener, 1936 1 10 0 0 
O. (Melanosmia) phenax Michener, 1936 0 4 0 0 
O. (Melanosmia) prunorum Cockerell, 1897 7 165 10 9 
O. (Melanosmia) unca Michener, 1937 0 8 15 0 
Osmia sp. 2. 1 9 0 0 

Mellitidae 
Dasypodinae 
Dasypodaini 

Hesperapis (Amblyapis) larreae Cockerell, 1907 0 32 0 0 
Hesperapis (Disparapis) sp. 1 27 40 0 0 
Hesperapis (Panurgomia) fuschi (Viereck, 1909) 1 144 1 0 


“core area’) made from 2000 to 2008 (Fig. 1). Not included were the samples made by 
pan traps at permanent sites in 2001. In this dataset, samples were included regardless 
of sampling technique (pan trap, aerial netting at flowering plants) or duration, and 
included collections made at the permanent sites and intervening sites (Tables 2, 3). 

We tested if results from the 2001 intensive pan trap study differed from the results 
from the same “core area” sampled from 2000-2008 by comparing the accumulation 
curves from both datasets after they were standardized for equal sample size and coverage 
following the approach by Chao et al. (2014). Predicted species richness was estimated 
after the two samples were extrapolated to 42,798 individuals, twice the abundance 
of the smaller sample (pan traps from 2001 at the permanent sites, Nj = 21,399 
individuals) as recommended by Chao et al. (2014) and implemented in the program 
iNext (http://chao.stat.nthu.edu.tw/wordpress/software_download/inext-online/). 
Confidence intervals and standard errors were calculated from 50 bootstrap replicates. 
Curves were significantly different if confidence intervals did not overlap. The 
contribution of vegetation types to bee species richness is not considered in this study. 
An earlier study found that there was little difference among vegetation types in bee 
species richness and composition, most likely because of the close proximity of the 
habitats and sites sampled (Minckley 2014). 
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Figure 2. Species estimates differ among single year and multi-year bee sampling. Comparison of sam- 
ple-sized-based rarefaction (solid lines) and extrapolation (dashed line) from the pan traps samples at 
permanent sites in 2001 and the samples made by netting and pan traps throughout the same area en- 
compassed by the permanent sites over 9 years (the “core” area). Calculations here are for Hill numbers of 
order q=0; orders q = 1 and q = 2 are in Suppl. material 1: Figures S1-S3. Bee species richness at the base 
sample size of 42,798 individuals (i.e., double the smaller reference size) is 406.0 for samples from the core 
area and 321.8 from the 2001 pan traps. Although bees were sampled from the same area, differences in 


sampling and years sampled gave different species richness estimates. 


Species richness 


To examine if the samples from the valley floor approximated the species richness of 
the fauna in a larger regional area, we compared species composition and estimated 
species richness to our complete dataset of bees from the San Bernardino Valley that 
includes all bee samples including two from nearby areas. One area extended 18 km 
east from the valley floor area up to desert grassland habitat. We refer to the collections 
from this area as from “peripheral” sites (Fig. 1). These samples were made during the 
months when samples were made in the valley floor and used pan traps and aerial net- 
ting. A second set of samples was made from 2009-2014, as part of a study of cactus 
reproductive biology by W. Radke, at a limestone outcrop dominated by desert vegeta- 
tion approximately 10.5 km north-northeast of the valley floor area (Fig. 1). Most, but 
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Figure 3. Estimates from all sampling regimes predict approximately 500 bee species in the San Ber- 
nardino Valley. Comparison of sample-sized-based rarefaction (solid lines) and extrapolation (dashed line) 
from samples made in the core area, core area combined with the 2001 pan trap sampling at permanent 
sites, and all samples made in the San Bernardino Valley. Numbers in parentheses are the number of in- 
dividuals in each dataset and the number of species. Calculations here are for Hill numbers of order q=0; 
order q = 2 is in Suppl. material 1: Figure S4. Estimated bee species richness at the base sample size of 
94,240 individuals (i.e., double the smaller reference size of the smallest dataset) is 499.3 species, 508.4 
and 516.1 species for the Core, Core + 2001 pan trap and All datasets respectively. 


Table 2. Sampling effort for bees among four areas in the San Bernardino Valley. Pan trap sampling fol- 
lowed the protocol of LeBuhn et al. (2007). 


Sampling years Sampling months Collections (date x site) Floral host sp. 
pans net 
2001 pan 1 Apr-Sep 848 848 n.a. na. 
Core area 9 Apr—Sep 1370 788 590 128 
Perimeter 9 Apr-Sep 131 109 22 61 
Limestone 5 Mar—Apr 44 29 15 2 


not all, of these samples were made in the spring months and included sampling by 
pan traps and nets. Tables 2, 3 show details on the number of collections and sampling 
protocols used for these two areas. 

Here, and as described above, the standardized curves of species interpolation and 
extrapolation and sample coverage were generated by the program iNEXT (http:// 
chao.stat.nthu.edu.tw/wordpress/software_download/inext-online/), and comparisons 
of datasets and estimates of species richness were done as recommended by Chao et al. 
(2014). Predicted species richness from samples made in the core area, core area plus 
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Table 3. Specimen collections of bees from four areas in the San Bernardino Valley. 


# individuals # species # unique species 
2001 pan 21403 277 23 (8.3%) 
Core area 47120 447 132 (29.5%) 
Perimeter 3253 221 12 (5.4%) 
Limestone 2375 77 12 (15.6%) 


2001 pan traps, and samples from all sites made in the region was estimated after the 
samples were extrapolated to 92,240 individuals, or twice the abundance of the small- 
est sample (core area collections, N oheerved = 47,120 individuals). Confidence intervals 


and standard errors were calculated from 50 bootstrap replicates, and accumulation 
curves were considered significantly different if confidence intervals did not overlap. 


Comparison to other studies 


Meiners et al. (2019) recognized that comparisons among long-term bee studies would 
be more informative if differences in the area sampled were considered. To standardize 
species-area relationships, they generated the species-area curve from 22 long-term bee 
studies in the United States and then ranked each study based on the percent difference 
between observed and predicted species density. Meiners et al. (2019) included the San 
Bernardino Valley study in their analysis but calculated species density using the area of 
the San Bernardino National Wildlife Refuge (1088 km”), instead of the area among 
permanent sampling sites (16 km’), and an earlier estimate of bee species richness 
(N = 383) reported in Minckley (2008) than known at present (N = 473). We used the 
same data and study sites included in the Meiners et al. (2019) study to compare new 
results on species density from the San Bernardino Valley to the other studies. Table 4 
in Meiners et al. (2019) with changes updated to reflect the findings in this study are 
in Suppl. materials 1: Table S2. 


Results 


Samples limited only to pan traps at permanent sites in 2001 yielded a total of 21,399 
specimens representing 277 bee species (Table 1) from 848 samples (day x site). Sam- 
ples over nine years and more collections (N = 1370) with both pan traps and netting 
at flowers (Table 2), and from the same area defined by the permanent sites, yielded 
more than twice the specimens and 447 species (Table 3), including 58 species not 
collected in the 2001 pan trap dataset. Of the total specimens represented when the 
two datasets are combined, 31% were from the intensive pan trap sampling in 2001 
and 69% were from sampling using nets and pan traps across multiple years. ‘The spe- 
cies accumulation curves from these two sets of samples differ significantly after 200 
individuals, as judged by non-overlapping confidence intervals (Fig. 2), despite having 
been collected from the same area (Fig. 1) and representing robust samples based on 
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estimators that include sample coverage (Suppl. material 1: Figures S1—S4). The core 
area samples rarefied to the number of individuals captured in the 2001 pan trap sam- 
ples (N = 42,798) yielded significantly more species (N = 437.6) than the 321.8 species 
predicted from the 2001 pan trap samples (Fig. 2). Overall, the observed species rich- 
ness from the intensive 2001 pan trap sampling yielded 62% of the species collected 
from the core area over all 9 years. 

Far fewer collections of bees were made from sites outside the core area and these 
amounted to fewer individuals and species of bees (Tables 2, 3). Nevertheless, the 221 
species from the peripheral and limestone samples are a large fraction (79%) of the 
277 species captured by pan traps in 2001 and about half the number of species repre- 
sented in the core area dataset (Table 3). Twelve species were unique to the peripheral 
site samples. The limestone outcrop samples were the smallest in terms of specimens 
collected and species (Tables 2, 3), however, 12 of the species were unique even though 
sampling was focused primarily on pan traps and bees were recorded from only two 
floral hosts (Table 2). 

Including all samples at all sites in the area defined by the permanent sites (2001 
pan traps and core area) yielded 473 bee species (Table 1). With the additional sam- 
pling made at the peripheral sites and the limestone outcrop there are 497 bee species 
known from the San Bernardino Valley (Table 1). Overall, the two datasets from the 
area defined by the permanent sites (2001 pan traps and core area) captured most of the 
bee species richness in the desert areas of the San Bernardino Valley; samples for sites 
peripheral to the core area increased the total number of bee species by only 12 species. 

Pooled across all sites in the region the collections total to 74,147 individuals and 
497 bee species (Table 1). Although the core plus 2001 pan trap dataset has many 
more individuals (N = 68,519) than the core dataset (N = 42,120), it results in a spe- 
cies accumulation curve that is significantly lower than the curve based on all datasets 
combined (Fig. 3). Up to approximately 42,000 individuals, the species accumulation 
curve from the core dataset is significantly lower than the curve based on all datasets 
combined. However, at greater numbers of individuals the curves are statistically in- 
distinguishable based on overlapping confidence intervals (Fig. 3). Extrapolations of 
species richness to 94,240 individuals predicted 508.4 species for the core area dataset 
and 498.3 species for the 2001 pan trap combined with the core area dataset and 
closely approximate the estimate of 516.1 bee species based on all samples made in the 
San Bernardino Valley (Fig. 3). The increasingly flattened accumulation curves as more 
individuals are sampled indicate these estimates of bee species richness are unusually 
robust (see Williams et al. 2001), and biologically reasonable, given that they closely 
approximate the number of bee species (N = 497) from all samples combined made by 
us in the San Bernardino Valley. 

Among the 22 long-term bee studies that Meiners et al. (2019) listed, three (Yo- 
semite National Park, Clark County, Nevada and Grand Staircase Escalante National 
Monument) reported more bee species and only one (Black Rock Forest Reserve, New 
York) was from a smaller area than the core area sampled in the San Bernardino Valley 
(Suppl. material 1: Table S2). The updated information from the San Bernardino Valley 
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Figure 4. Higher bee species density in the San Bernardino Valley than that reported from other long- 
term studies. Species density of bees for large-scale studies in North America determined by species-rich- 
ness relationships following Meiners et al (2019) with updated information from this study. a) the trend- 
line indicates expected change in species richness with area sampled based on 22 studies (listed in Table 4 
of Meiners et al [2019]) plus the San Bernardino Valley (black dot circled in red). b) is the difference in 
bee species observed relative to the prediction from the trend line in a). Updated information from the 


San Bernardino Valley alters the rank of species density from that shown in Fig. 5 of Meiners et al (2019). 


changed the species area curve and the rank order of species density for the 22 long- 
term bee studies as judged by the percentage change between observed and predicted 
species per area (Fig. 4a, b). Meiners et al (2019, see their fig. 5a, b) concluded that the 
six sites with the greatest percentage increase over that predicted from the species-area 
curve was Pinnacles National Park, California followed in descending order by Grand 
Staircase Escalante National Monument, Utah, Yosemite National Park, California, 
Carlsbad Caverns National Park, New Mexico, Clark County, Nevada and the San 
Bernardino Valley, Arizona/ Sonora. In our study with the new information reported 
herein, the San Bernardino Valley, Arizona/ Sonora is the highest ranked study fol- 
lowed in descending order by Grand Staircase Escalante National Monument, Utah, 
Pinnacles National Park, California, Yosemite National Park, California, Clark County, 
Nevada, and Carlsbad Caverns National Park, New Mexico. With these changes, bee 
species density from the valley floor samples in the San Bernardino Valley is greater 
than any other area in North America surveyed to date (Fig. 4b). 


Discussion 


The San Bernardino Valley of the Chihuahuan Desert is extremely rich for bees. An 
area of 16 km’ that varies in elevation by less than 120 m yielded 277 species in 
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one year of sampling and 473 species over nine years. The four bee surveys in North 
America that have reported more bee species than are known from the San Bernardino 
Valley (N = 656 species, Great Staircase Escalante National Monument Utah; N = 598 
species, Clark County Nevada; N= 554 species, Yosemite National Park California; N 
= 450 species, Pinnacles National Monument California) are much greater in area and 
topographic complexity (Suppl. material 1: Table S2, Meiners et al. 2019). When the 
effect of area alone is standardized, the density of bee species in the San Bernardino 
Valley, is the highest reported from any long-term study of bees in North America. 

For questions that rely on accurate estimates of species composition and richness, 
this study highlights the importance of using a variety of sampling protocols and long- 
term sampling. Our dataset based on samples from one passive sampling protocol (pan 
traps) in 2001 that recovered more than 21,000 specimens yielded 93 fewer species 
than we found from samples made over nine years that included pan traps and aerial 
netting at flowers (Table 3). If we assume that the collections from 2001 pan traps and 
core area are complete measures of species richness, pan traps underestimated species 
richness by 38%. The significantly lower species accumulation curve from the 2001 
pan trap samples than the long-term dataset from the same area (Fig. 2) is further evi- 
dence that surveys based only on pan traps underestimate bee species richness. Biases 
of pan trap sampling for bees have been well documented (reviewed in Packer and 
Darla-West 2021) but not from a long-term study in a bee diversity hotspot. 

There are numbers of reasons why species richness estimates for desert bees are 
sensitive to sampling duration that broadly apply to biodiversity studies of most de- 
sert insects and bees in general. For one, many species, especially those bees that are 
specialists of one or several host plants, emerge as adults for a brief period per year, 
are small-bodied with limited dispersal capabilities and occur patchily (Minckley et al. 
1999; Kazenel et al. 2020). Many species are rare: singleton (represented by one speci- 
men) and doubleton (represented by two specimens) species represent 83 (16.7%) and 
52 (10.5%) of all species in our full dataset, respectively. Finally, bee populations vary 
dramatically between years in response to changes in floral resources (Crone 2013) and 
propensity to remain in diapause (Minckley et al. 2013). Variation of species abun- 
dance taken at the same location across years is often marked (Stubblefield et al. 1993; 
Meiners et al. 2019). An extreme case occurred in spring 2006 when a drought resulted 
in the complete absence of 11 of 46 bee species that had been abundant in the San 
Bernardino Valley in the near-normal rainfall years before and after (Minckley et al. 
2013). Of the absent species, seven were pollen specialists of plants that require rainfall 
to trigger growth and flower (annual species and creosotebush [Larrea tridentata]), 
and four were pollen generalists. In contrast, the five pollen specialist bee species of 
mesquite (Prosopis velutina), a deep-rooted perennial that blooms regardless of winter 
rainfall, were active during the drought. 

The ability to undergo facultative diapause is common among desert bee species 
(Hurd, 1957; Rust et al. 1988; Danforth et al. 2019) and suggests a bet-hedging life 
history strategy is a common adaptation for persisting in these areas (Minckley et al. 
2013). Danforth (1999) showed in a laboratory experiment that some Macrotera portalis 
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(Timberlake, 1954), a specialist of mallows (Sphaeralcea spp), emerged in response to 
moisture in the late-summer (when it is active) but others remained in diapause for one 
or two years. Desert plant species vary tremendously in their propensity to emerge and 
bloom (Bowers and Dimmitt 1994; Bowers 2005), even in common gardens (Gremer 
and Venable 2014; Gremer et al. 2016), which suggests that activity of the specialist bees 
that depend on different hosts will differ in their activity between years and across habitats. 
Sampling that is only occasional through the season or widely-spaced across an area may 
not detect many resident bee species. In the San Bernardino Valley, there is 80% turnover 
of pollen specialist species every two weeks in the spring (Minckley 2008), suggesting that 
samples for these and other short-lived solitary bees should be more than bimonthly. The 
samples in this study made over 9 years and using pan traps and nets minimize these biases 
and account for the higher estimates of bee species richness than from very intensive pan 
trapping in one year even when both datasets are from the same location. 

The species richness estimates using extrapolation methods from the samples on 
the valley floor suggest approximately 500 bee species occur in the San Bernardino Val- 
ley, of which we captured 497 when all samples in the region were combined (Fig. 3, 
Table 1). Extrapolations of all samples suggests species richness is 516.1, not substan- 
tially more than the number of species observed, and approximately 14% of the 3562 
species known for all of North America north of the United States-Mexico border 
(Ascher and Pickering 2020). The number of bee species from the San Bernardino 
Valley and their relative abundance sets a reasonable expectation for communities in 
other undisturbed desertscrub/ desert grassland habitats of the Chihuahuan Desert 
and can be usefully compared to areas with more severe human impacts. For example, 
a compilation limited to described bee species from southeast Arizona in the American 
Museum entomology collection and from the San Bernardino Valley yielded 540 bee 
species when desert scrub, mid-elevation grassland and high montane ecosystems were 
included (Minckley and Ascher 2013). Considered individually, the three desert areas 
north of the San Bernardino Valley had between 150-180 species but did not include 
some genera that had yet to be curated in the American Museum of Natural History 
collection. Thus, further taxonomic and curatorial work should substantially increase 
the bee species richness from these well-collected sites. Similarly, a 12-year study of 
bees on the extreme northern edge of the Chihuahuan Desert in New Mexico (approx- 
imately 510 km NNW of the San Bernardino Valley) captured 308 bee species using 
a modified pan trap design (Kazenel et al. 2020). Because bee species richness is often 
underestimated by pan traps, the reported number of species suggest this bee fauna is 
much richer and may be comparable to that in the San Bernardino Valley. 

Museum records, taxonomic revisions, and surveys have long suggested that bee 
diversity runs counter to that of the plants they pollinate in that bees have peak diver- 
sity in the western hemisphere in the Chihuahuan, Mojave and Sonoran deserts and 
in the eastern hemisphere around the Mediterranean Sea (Michener 1979, 2007). The 
repeat sampling from a defined area in the Chihuahuan Desert provides empirical data 
that supports this counterintuitive biogeographical pattern. However, more sampling 
of bees is needed in tropical and subtropical areas, as well as a more complete picture 
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of herbaceous plant biogeography. Plant biogeography still rests largely on what we 
know about tree species, and when short-lived plant species are better documented 
(Massante et al. 2019) their distribution may more closely mirror the distribution of 
their primary pollinators, the bees. 
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Tables $1, S2, Figures $1-S4 

Authors: Robert L. Minckley, William R. Radke 

Data type: tables and figures 

Explanation note: Table $1. List of identification keys used for each of the genera 
sampled in the San Bernardino Valley. Table $2. Table 4 from Meiners et al (2019) 
of bee species density for native bee inventories with at least 100 species in natural 
or semi-natural areas across the United States. Information here is updated for the 
San Bernardino Valley (this study). Figure $1. Comparison of sample-sized-based 
rarefaction (solid lines) and extrapolation (dashed line) from the pan trap samples 
at permanent sites in 2001 and the samples made by netting and pan traps through- 
out the same area encompassed by the permanent sites over 9 years (the “core” area) 
following the approach of Chao et al (34). Bee species richness at the base sample 
size of 42,798 individuals (i.e., double the smaller reference size) is 406.0 for sam- 
ples from the core area and 321.8 from the 2001 pan traps. The two top graphs 
are sample-size-based rarefaction and extrapolation curves for the pan trap samples 
from 2001 and the samples made in the core area from 2000-2008. The bottom 
three graphs are comparisons of samples from the pan trap samples from 2001 and 
the core area for each order of q. Although bees were sampled from the same area, 
differences in sampling and years sampled gave different species richness estimates 
for each order of q. The first graph in lower row is the same as Fig. 1 in the paper 
and presented here for completeness. Figure $2. Plot of sample coverage for rarified 
samples (solid line) and extrapolated samples (dashed line) as a function of sample 
size for the pan trap samples at permanent sites in 2001 and the samples made by 
netting and pan traps throughout the same area encompassed by the permanent 
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sites over 9 years (the “core” area). 95% confidence intervals calculated from 50 
bootstrap replications are not prominent because of their small size. Both sam- 
ples have greater that 95% coverage and broadly overlap up to their reference sizes 
(21, 403 individuals for 2001 pan trap samples, 47,120 for core area samples) and 
when extrapolated to twice the reference sample size. Figure $3. Comparison of the 
coverage-based rarefaction (solid line) and extrapolation (dashed line), for the pan 
trap samples at permanent sites in 2001 and the samples made by netting and pan 
traps throughout the same area encompassed by the permanent sites over 9 years 
(the “core” area) for all three orders of q. In each of the three graphs the lower line 
is the 2001 pan trap samples and the upper line is the samples from the core area. 
Note that the lines do not intersect at any sample coverage. Figure $4. Comparison 
of the coverage-based rarefaction (solid line) and extrapolation (dashed line), for 
the samples made in the core area (middle line), core area combined with the pan 
trap sampling at permanent sites (blue line on bottom), and all samples made in 
the San Bernardino Valley (top line) at q =0. The confidence intervals are very small 
they do not show on the graph. Sample coverage is above 95% for diversity values 
below 200 species for the datasets and the extrapolations overlap broadly at higher 
sampling coverage indicating all three datasets predict very similar bee species rich- 
ness in the San Bernardino Valley. Note that the analyses at other orders of q are not 
shown because the results match that shown here. 

Copyright notice: This dataset is made available under the Open Database License 
(http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License 
(ODDbL) is a license agreement intended to allow users to freely share, modify, and 
use this Dataset while maintaining this same freedom for others, provided that the 
original source and author(s) are credited. 
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